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ABSTRA.QT 

The modern c aitrol theoretic approach to the 
design of the Voltage Regulator and Speed Governor for 
a Synchronous Generator is given. An accurate linearized 
state space representation of the synchronous generator 
is formulated for ?/hich an optimal output regulator is 
designed. To implement this control a dynamical 
compatible observer which reconstructs the state vector 
is employed. The superiority in the performance of the 
system with optimal output feedback regulator compared 
to that with a conventional one is conclusively 
established. 
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OHAPTEE I 


IHTRODUOTIOF 


In recent years there has been a great deal of 

interest in the raathenatical modeling of synchronous 

generrtors along wath their c cut rolling equipment mainly 

from the stand point of dynamic stability analysis, and 

the study of system performance during disturbance or 

abnormal c onditions°’ , more recently there have also 

boon attempts at the application of optimal control 

8 

theory C'Oncepts to the design of the regulators , 

In this thesis, a new approach to the design of 
the volta-j^e regulator and the speed governor for the 
synchronous generator is given-. This approach is diffe- 
rent from the conventional design procedure, wherein the 
designer chooses a prior'^configuration for his controller 
and tries to adjust the various parameters to meet his 
requirements in s o far as possible. In this approach an 
optimal feedback control law is first obtained for the 
synchronous generator which is then implemented by means of 
a dynamical compatible observer, the optiwial control 
la.y^ is cascaded ■■ with the observer configuration resul- 
ting in the desired Optimal Output Regulator. 
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Oiaapter II deals with the mathematical represen- 
tation of the synchronous generator in a form suitable for 
application of the techniques of modem control theory. 

The linearized state space model cf the system developed, 
forms the basis for the investigations in the later 
chapters. 

The conventional design procedure of the regulators 
is discussed in Oliapter III and the performance of a genera- 
ting system with a typical configuration of Voltage Regulator 
and Speed Governor is investigated. 

In Chapter IV, an infinite time linear output 
regulator problem is formulated.' The choice of this parti- 
cular criterion is mainly dictated by the ease in the 
implementation of the resulting optimal controller. The 
optimal control law is obtained as the solution of the matrix 

Riccati equation, which has been solved by the technique 

1 6 

of Successive Approximation , ' f or different operating 
conditions and perfoimaice indices- The implementation of 
the optimal control law requires the entire state vector 
for the purpose of feedback, a condition which is difficult 
to satisfy in many practical situations. 

To obviate this difficulty an observer system".- 


can be employed, the design of which is treated in 
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Chapter Y, it is shown that the observer system reconstructs 
the state vector from the available system inputs and outputs 
and the reconstructed state vector is used in the implementa- 
tion of the optimal feedback control law. The perfoimance 
of the Optirfial Output feedback regulator is investigated 
and compar'ed to that of the conventional regulator. The 
practical realisation of the former is also indicated. 



CHAPTER II 


STATE SPACE MODEL OP THE SYNOHROROUS GENERATOR 

2 , 1 Introduction ; 

1 -5 

In the earlier studies , the mathematical 

rej,jresentation of synchronous, generator was based on the 

generalized theory of Kron or on the operational impedance 

4 5 

aj,.proach of Park and Concordia . Recently, inspired by 

6 

the developments in control system theory, Laughton , 

7 8 

TJndrill and Yu etal have extended the study of dynamical 
system behaviour, from the analysis of differential equation 
properties as exhibited by single higher order equations to 
the more general vector space methods. 

In this chapter a state space model of the 

synchz'onous 'generator is formulated. The set of state 

variables chosen to describe the synchronous generator are 

the flux linkages of the direct and the quadrature axis 

armature and rotor circuits. These set of variables offer 

a definite advantage over the set of current and voltage 

8 

variables chosen by Yu etal . Since the flux linkages in a 
closed circuit cannot change suddenly, hence the initial 
conditions for any system disturbance can be found easily. 
The complexity in carrying out the analysis on the non- 
linear representation of the system is simplified by 
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linearization of the system differential equations about an 
operating point. 

2 . 2 Hathematical Representation of Machine and 
Transmission Line ; 

A generating system, where the generator feeds a 
bus of known voltage and phase angle via a long transmission 
line with an addit'ional shunt load at its sending end, as 

7 

studied by Undrill , is chosen for investigation (see Fig. 2.1). 

A complete description of the'dynamio behaviour of 
the synchronous generator requires consideration of its 
electrical and mechanical characteristics. The necessary 
mathematical statements as given by Park^ in per unit and 
extensively used in literature are summarized in the following 
paragraphs. The assumptions made are; 

1 . The stator windings are sinusoidally distributed 
around the air gap as far as mutual effects between them and 
the rotor are concerned. 

2. The stator winding self and mutual inductances vary 
sinusoidally as th€.TOtor moves, and are of the form 
a+bCos2©5 a-j+bCos2(0 + ^) , respectively, where a, and b 
are constants. 

3. Hysteresis effects are negligible. 

4. Electrical transients on the transmission system 


are neglected. 
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5. line parameters do not change with instantaneous 
generator frequency. 

■ 6. Leakage reactances of all the windings are independent 
of the saturation state of the iron. 

7. The leakage fluxes do not contribute to the iron 
saturation, which is therefore determined by the mutual flux. 

8. Mutual reactances between the direct axis armature 
winding and the two direct axis rotor circuits are made equal 
by S'uitable choice of the base rotor currents, and owing to 
the proximity of the amortisseur winding to the air gap, it is 
assumed that the mutual reactance betw^een the two direct axis 
rotor circuits is equal to that between these circuits and 
arm ature . 

9. Only those modes of operation that do not require 
zero axis variables are considered. 

E lectrio.al Equations ; 

A typical layout of the rotor field winding and 
stator armature', windings of a 3~phase, 2 pole synchronous 
generator is shown in 5'ig.2.3. The path of the current flow 
in the rotor is represented by windings in direct and 
quadrature axis. The voltage and current equations of the 
synchronous generator are written for the individual circuits 
of the generator and mutual effect between them. These 
equations contain terms which are functions of rotor angle 0. 
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The difficulty in solution is overcome by Park's transformation. 
It resolves the statox’"' quantities into components along the 
direct and quadrature axes. The resulting equations obtained 
for a balsaaced J^phase performance are; 


Di x;pct axis flux linkage ; 

Pu ~ ^ffd^fd + ^ad^kd " ^ad^d 

f-'d ■ ^ad^fd ^ad^kd " 

{^’’kd "" ^ad^fd ^kkd^kd “ ^ad^d 

Quadrature axis flux linkage ; 

f q ~ ^aq^kq “ ^q^q 
fkq ^kkq^kq “ ^aq^q 
D irect axis voltage - 

"'^fd^ ^ ^ fid ^fd^fd 

^d= 

0 = — 

Quadrature axis voltage ; 

®q ~ ^ P% *■ ^^q 

^ ^ ^ PKq ^kq^kq 


( 2 . 2 . 1 ) 

( 2 . 2 . 2 ) 

( 2 . 2 . 3 ) 


( 2 . 2 . 4 ) 

( 2 . 2 . 5 ) 

( 2 . 2 . 6 ) 

( 2 . 2 . 7 ) 

( 2 . 2 . 8 ) 

( 2 . 2 . 9 ) 

( 2 . 2 . 10 ) 
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Additional Equations: 

B lcctrioal torque in air gap ; 

Reactive power at terminals ; 

Q = eqi^ -- e^^i^ 

G-enerator terminal voltage ; 

2 2 2 

®t = ®d ®q 

T'echanioal eqaation ; 

"i = + K^(pS) +At 

Terminal constraint equations; 


( 2 . 2 . 11 ) 


( 2 . 2 . 12 ) 


(2,2.15) 


(2.2.14) 


The nodal equation for a T-equivalent of the transmission 
system, shown in Eig.2.2, is 


e - e’ 




e'Y^ + 



V 


or 1 (Z^+Z^.Z^Y 3 .Z^) = e(YgZ^+1)-T 


<2.2.15) 


Considering the d and q axes of the machine as the reference 
axes, the phase representation of terminal current and voltage 
is given by 

i = ijj 4- jiq 

e = e^j + jSq 
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Substituting the values 

(i^+jiq)((Rg+jXg)+(R^+jX^)(Rg+jXg)(Gg+jBg)+(R^+dX^)) 

= (e(j+jeq_)((Rg+DXg)(Gg+jBQ) + 1 )-(Vjj+jVq) 

Equating real and imaginary terms 

(2.2.16) 

|j ^d*'" ~ ^ ^d "^^q ^q (2.2.17) 

y,here - (X^B^R^+E^B^X^+X^G^X^) 

r = X^ + + EgB^E^ + X^H^G^ + EgG^X^ - X^X^B^- 

C< = 1 + E^G^ - X^B^ 

P = E„B. + X„G„ 

I c t; ye 

Rearranging the equations (2.2.16) and (2.2.17) 


I) 

CD 

ei^ + 

fiq + 


+ hVq 


(2.2, 

.18) 

®q = 

-fi^ + 

eiq 

- 

+ 

g^q 


(2.2, 

.19) 

where 

e = 

~-T“ 

y- p> 

9 

f 

_ H 
~ ? 

- cca^ 

£^2 




OC + j 




oc 

+ p 



g = 

<x 

2 


9 

h 

^ ~ 

— ^ 




z<. + }■ 

6 



oC 

+ S 




To apply the techniques of modern ccntrol theory, if 
is convenient to represent the linear constant co-effioient 
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system In' the form x = A x + B u , and the output equation 
?v’ritten as Z = 0. vector u comprises of the input 

variables to the system. The variables of vector u comprises 7 

of the input variables to the system. The variables of 
vect-or Z. available outputs of the system, which is 

being controlled. The vector x denotes the state of the 
system. The above set of equations (2.2.1) - (2.2.19) are 
manipulated to obtain a state space representation, the states 
being 5’^ .^2 and "tlie output 

vari ables . 


Substituting the values of the currents in terms of the flux 
linkages in equations (2.2.1) ^•'•(2.2.19) results ,in 


Pffd 


= Wq( 


- ffd^ 


X. 


f1 


■^fd ®fd 

aa 


(2.2.20) 


P'% 


= W 


o^®a- ^ 


W 


- fa) 


'a1 


) 


= »o( i" + 




X 


a1 


) 




w 




‘•kql 


pfA = - 

pd^ = 62 


w. 


^kd ^Tad ~)^kd^ 

^kdl 


( 2 . 2 . 21 ) 


( 2 . 2 . 22 ) 


(2.2.23) 


(2.2.24) 


(2.2.25) 
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pS 


-^2 




= - Tg - KjSg -at) (2.2.26) 

^ ,vSlnS,.hvCo.S, 

^a1 ^a1 ^ ' 

. . .(2.2.27) 

^a1 ^a1 ^ ^ 

. . .(2.2.28) 


'a1 


1 / ^d ' ^^fd ^kd v 

^1 ^a1 ^f1 ^a1 


‘)ir' 

i aq. 


= i_(i^ + 

^2 ^kq1 ^a1 


ciinCi 


K, = 




1 


X 


1 , I 

+ t"! + 


ad 

1 


^f1 

1 


X 


■k1 

1 


X 


a1 


X 


aq 


^kdl ^a1 


Substitution "f'ad ’ ®d equations 

(2.2.20) - (2.2.28) results in 




7fd 


^o^fd^ 1 


x^1 ^K^x^., 


- 


fd o r^/f- , fd 0 r,i/' - 

Tfi + =;^rkd 








Pfa 




I’vJv+e) 


WQ(r+e) 

^l^fl^al^^^ ^ ^a1 ' '^l^al 


-(■ 


- uK 


(2.2.29) 
w_(r+e) 




'P\d 


+ (w + ^( 


' - U)F„ . 1V.r-r. 


^a1 Val 




+ 'WQgvSinO.j + w^hv0os5'^ 


( 2 . 2 . 30 ) 
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pV; . = . ^o^kd o/. Vkd. 1 


^^0 . w f 

KiXgiX^^7fd ■*■ + x77^^ ~ f 


. . .(2.2.31) 




^al^kd l 


w (r+e) 


'a1 '^2 




<3. ^2^a1^kq1 ' 


- v;?^hvSin5^ + w^gvGosd^ 


(2.2.32) 


w p, 
o-^kq 




2^k.q1^a1 ^ ^ ^kql 


(2.2.33) 


(2.2.34) 


rdi - K,s, -'^., . titM ti'^i 

rn ^ rr- ? 

2 al ^a^al^fcll *^1^a1 


+ ka?fd , f q , 


(2.2.35) 


-22 
®d + ®q 


(2.2.36) 


2.3 Linearized l^/rndpi ^ 


(2.2.37) 

(2.2.38) 


The system under consideration is basically a 
Jnlinear system. Its performance is investigated by means 
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of linearized model. This form of analysis has been 
extensively used in literature . It retains the first term 
of the Taylor series expansion of the equations about the 
operating point. The equations (2.2.29) - (2.2.38) are 
written for perturbed values of voltages and fluxes (such as 
‘ Then subtraction of the original equations and 
neglecting of second order terms, yields the linearized Model. 
The details of this'are given in Appendix (a)« 


Waen linearized the system eOrations can be written 
in the following standard form 


1 


1 d 

i 


— 

®ii 

^12 

®13 

0 

0 

0 

0 


m 

^1 

0 

i pA|/-d 

j 


®21 

®22 

®23 

^24 

^25 

®26 

^27 

llfi 


0 

0 

P-Af kd 


®31 

®32 

®33 

0 

0 

0 

0 



0 

0 


= 

^41 

^42 

®43 

^46 

^45 

^46 

®47 

Afi 


0 

0 

1 1 

pAfi., 

! 

0 

0 

0 

®54 

^55 

0 

0 

hfc 


0 

0 

1 


0 

0 

0 

0 

0 

0 

1 

AS, 

' i 


0 

0 

pA52 


®71 

^72 

a? 3 

®74 

®75 

0 

®7T 

AS 2 ! 


0 

^72 


"fd 


At. 


(2.3.1) 


where/Vj^^ ,l\% are the state 

variables of the system? andAE^^,AT^ are the forcing inputs 
to the system. 
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Tile output equations of the system are 


|Aet 


°1 1 

C 12 

°13 

3=14 

^15'^16 

0 

h5i 

= 

0 

0 

0 

0 

0 

1 

0 

AS 2 


0 

0 

0 

0 

0 

0 

1 


ATfd 

Afd ! (2.3.2) 

Afkd 

A% 

Afkq 

A S 1 
AS, 


■,a= 6 re Ao^.,A 5 ,,A 62 are the aTailable outputs of the system. 

• 1 ) 


a 


a 


"'^o^fd , 1 


11 


^f1 Kl^fl 


''^o^fd 


^ ^o^fd 
‘12 K.|X^^x^., 


13 


w^Cr+e) 

'21 = K^Xf^x^^ 


'22 


Wo(r+e) ^ w^(r+e) 


^a 1 


25 K^x^.,Xj^^., 


fWo ^ fw^ 

3o/1 (xr~tr ~ 


■24 - "0 


25 - 


®26 = e«o 


tGosS.| Q-hw^TSin^ q ®27 'If^c 


qo 



15 


a 


= 

Kl^kdl'^fl 






®33 “ X 


^"o^kd^ 1 


kdl ^l^kdl 


- 1) 


fw. 


3/1-1 — ”■ 


"41 


^l^al^f 1 


®42 ~ 


w„f . 

0 / ^ I 

- w^ -f -r — (1- 


0 X 


a1 


^l^al 


■) ® A ■3; 


fw. 


45 


a 


44 


w (r-^e) ^ 

-(r4- -1 ) 


"al 


■^2^31 


O '/} IT — 


"45 


w^Cr-f-e) 

^2^a1^kq1 


®46 "" “ 


10 s‘''o’ 


ho ®47 


= -fe 


do 


a = 

54 ICoXi,«mX, 


"2^kgl^a1 


a 


55 


-(xT^ “f) 


^kql '^2^kq1 


a , - 

^71 “ T K 


m'^'l^al^f 1 


aj2 - 


1 / ^o ykqo _ l^qo > 

'^m^al ^2^a1 '^^2^kq1 ^l^al 




"^5 '^m^l^al^kdl 


a 


74 


"I /^do Tdd l^do "fedo 


T x^T^k^^x . 
m a1 2 a1 


^l^al ^l^fl ^l^kdl 


) 


^75 


'^m^2^a1^kq1 


377 * 


T„ 


m 


*>11 = 


"o^fd 


X 


ad 


n2 ^ T . 


m 


= ^°<lo - 

•11 


•12 


^o^al ^l^al 


1)(eedo-feqo) 



16 


Cl 3 


- -"go 


'14 
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Cl 5 


®o^2^a1^k<l1 


=16 = -4;('"'=°=^10(S®do-*’®qo^ 
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2 .4 System Per-formanoe ; 

Por the purpose of analysis, a particular system 

7 

with constants as given by Undrill is taken and the relevant 
data are given below. 

Machine Data 


^ad 

= 1 .0 

= 1 .2 

^aq - 

0.6 

^ffd 

= m 

^kkq 

ii 

O 

CO 

^kkd ^ 

= 

0.8 

^a1 

= 0.2 

^f1 

= 0.1 

^kdl = . 

^kql " 

0.2 

^kd 

= 0.02 

^fd 

= 0.0011 

r = 0.01 

11 

0.04 

^m 

= 6 s 

Ki 

= 1 .0 

= 314.0 





ansmission Line Lata 





P^e 

= 0.05 

-•R^ = 0.0 

^6 = 

= 0.30 

^e = 

0.0 


= 0.0 

= 0.0 

6 






The operating point data calculated with an assumed 
saturation factor of unity are 

E^= 1 .0 = 0.5406 e^ = 0.8403 
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Tg = 0.6S41 6^ = 70. 00 = -0.6539 

Q = -0.5276 V = 1 .0 

These data are generally used in this thesis for 
the above mentioned operating conditions. The above system 
without an/ regulating eluipment is simulated on a digital 
computer and its response for a perturbation inAj^^j^O.Opp. u. 
is shown in ]j'ig.2-.4. Runge Ivutta fourth order numerical 
integration technique (see Appendix (B)) is used for a step 
size of 0.0025. The step size chosen is smaller than the 
srasllest time constant of the system. The response is rather 
poor, since even after the end of 5.0 seconds of real time,, 
the outputs A e^ and A 5^ continue to rise, although the 
output A ^2 tends to fall. The eigenvalues of this unregula- 
ted system are -0.046488, -0.26222, -0.63385, -27.826, 
-34.041, -71 .616 + 3635 . 50 . The negative real parts of the 
roots indicate that this system is stable. 

2.5 Conclusion : 

The state space model developed is an accurate 
representation of the system, since all the system parameters 
are included. This representation has an advantage over the 
other simplified models considered in literature. The 
simplified versions are usually , obtained by omitting the 
terms p^^ and p|/“ and sometimes the damper coils are assumed 



nonexistantj resulting in Ihe removal of rows and columns 
corresponding and'^j^q in the matrix representation. 

7 

The model is superior to that derived by Undrill , 
in the sense that no matrix inversion is involved, further 
it is ideally suited for the direct application of modern 
control techniques. The inherent poor response of the 
system obviously suggests the incorporation of a regulator 
in the system to improve the response. 
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CHAPTER III 

PSRPOMFGE A GOKVENTIOMAL VOLTAGE REGULATOR 

AND SPEED GOVERNOR 


3 . 1 introdaction ; 

Moaern generating sets are generally operated under 
less stable conditions than in the past. This combined with 
the current trend towards lower inertia constants and higher 
transient reactances renders the generators more likely to 
become unstable during fault conditions. For a fixed exci- 
tation, thex'e is a correspond irg constant value of the power 
limit of the machine. By the use of continuously acting 
regulating equipment, usually regulating the terminal voltage 

and speed to some fixed reference, the transient stability 

1 0 

limit of the generator can be improved 


In this chapter the conventional design procedure 
of regulating equipment is discussed. The perfoimance of a 
generating system with a typical configuration of Voltage 
Regulator and Speed Governor is investigated. 


•2 Conventional Design Procedure ; 


The conventional design procedure starts with the 

selection of a suitable configuration of the regulating and 

stabilizing equipment. Typical configurations are given by 
1 11 12 

Shackshaft , Meerov , Aldred etc. The system analysis is 


1 
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then carried out assuming the values for the time constants 
of the regulating and stabilizing equipment. The maximum 
overall gain of the loop is determined by the stability 
considerations, the most popular technique being the 
l-partition . The gain is then distributed to the various 
components in the control loop depending upon their physical 
feasibility. The system response is then observed for its 
transient response and steady state error. Depending upon 
the response, the gains and time constants are further 
adjusted, 

3 . 3 System Performaice v/ith Conventional Regulator : 

A fairly common representation of a continuously 

acting Voltage Regulator with derivative stabilizing 

transformer and a conventional speed governor is the one 

1 4 

used by Prabhashankar etal , shown in Dig. 3.1 •' nci 3*2- The 
control system equations are 

Voltage Regulator ; 

pEfd = " ^s^ - (3.3.1) 

Ph = - ''P (3.3.2) 

s 

■Speed Governor 

p(A ^ 1 ).= Alg 

p( A T 2 ) = - (Tg+Tj^)AT2 - At-, ) 

g h. o 


(3.3.4) 
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The above equations ( 3 * 3 .'i) -■ ( 3 - 3 . 4 ) are linearized 
and augumented to the system model, resulting in a state space 


representation of the form 
X = W X 

and X = I. 2 

where W is a 11 x 11 matrix and V is a 3 x 11 matrix, the state 
vector X is given by 

and the output vector x is given by 

y = ( Ae^A 6^ A62)^ 


The state space model of the system with the regulator is 
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The output equations of the system is 


p“‘®t| [®11 ®12 ®13 ®14 ®15 ®16 ® ^ ^ ^ ^ /iffd 


poi 


0 0 0 0 1 00000 


1^02 jO 0 0 0 0 01000 0 

‘ — I — 

where 

_ _ 1 f\ 

®18 - X.. ^710 ““T” , 

'^0. in AC*: 

il02 
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0 , S 

hh 
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The other elements are the same as defined in section 2.3 • 

For an initial perturnation in^'"^^ = 0.05 p.u. the 
.output response of the system with the above regulator is 
plotted in Fig. 3. 3 < The regulator parameters are 

Tj^ = 0.50 jUg = 0.04 Tg = 1.00 

Tq = 0.50 jj. = 5.00 jUg = 5.00 

These parameters are obtained by trial and error, so as to 
have a reasone.ble settling time and overshoot in the tran- 
sient response of the system, for an initial perturbation in 
the state variables. 

The oscillatory response of the system dies down 
in 10.0 seconds of the occurence of the initial perturbation 
. The oscillatory response is due to the dominant 
pair of complex eigenvalues of the system. The variation in 
parameters of the regulator has a marked effect on the 
response of the system. An improper choice of the gains and 
the time constants can render the system unstable. 

3.4 Oonolusion ; 

The conventional regulator design has a drawback 
in the arbitrary choice of the regulator configuration and 
the cut and try procedure involved in the selection of the 
values of the parameters to meet the specifications such as 
overshoot, decay time and stability limits. Even though 
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these specifications are achieved, it may not represent the 
best possible design. 

Through the choice of the appropriate performance 
index it becomes possible to impart the desired features to 
the 33/ stem transient response. Furthermore, it becomes 
possible to explicitly take into account any constraints 
involving the control and output variables of the system. 

In all eases the optimal control law derived exploits the 
maximum number of 'degrees of freedom' (.Horowitz) that are 
Eivsilable. In the next chapter the design of an optimal 
regulator for the power system under consideration is 
presented. 
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CHAPTER 11 

DESIGN OP Al'T OPTIFJM VOLTAGE AND SPEED REGULATOR 


4 . 1 latroductioi'i 

The design of a regulator for the synchronous 
generator, both for the control of the terminal voltage and 

speed, has been bestowed much attention in literature. 

8 

Recently, Yu etal have discussed an optimum output regulator 
for a synchrorjous generator system, already equipped with a 
voltage regulator and a speed governor. The output of the 
optj.mum controller is fed to two different summing junctions 
of the exciter voltage regulator and the speed governor 
through- a transfer function and a servomotor respectively. 

Thus, the problem considered by them is a "specific 
optimal" control problem, ■•/herein the structure of one part 
of the controller comprising of a conventional voltage 
regulator and a speed governor is assumed to be fixed and 
the optimal controller which forms the remaining part of 
the controller is obtained by treating the synchronous 
generator with the voltage regulator and speed governor 
as the equivalait plant to be controlled. The inclusion of 
conventional regulators is superfluous as the optimum 
controller is itself capable of providing the necessary 
control. This latter approach is adopted in this investi- 
gation. 
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The svstea is formulated as an infinite time 

output regulator control problem. The theory of the 

infinite time regulator problem for a linear time- invariant 

plant with a quadratic objective functional, with no explicit 

c onsjtraiiits involving the control and the output variables 
1 5 

is well known . for this class of problems, .the optimal 
controller obtained is linear and time-invariant and can be 
synthesised using the state vector feedback. 

for a linear regulator problem, the optirn?’! control 

law obtained requires the solution of a matrix Riccati 

equation. For higher order systems analytical solution to 

the resulting matrix Riccati equation becomes difficult. 

Therefore, the solution of the matrix Riccati equation is 

customarily obtained, by numerical techniques. The method 

1 6 

of Successive Approximation is fairly popular and is 
described in the following sections. 

4 . 2 Formulation of the Optimum Control Problem ; 

In this Section, the unregulated power system 

considered in Section 2.2 is formulated as an optimum 

1 5 

output regulator problem . The cost functional which is 

to be minimized is assumed to be of the form 

00 

J = i ^ R £ + H.) (4 .2 .1 ) 

where the matrix Q and the nonsingular matrix R are the 
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weightages associated vvith the output and the input 
variables and these are generally chosen to be symmetric 
matrices. The vi/eightages assigned -.'depends upon the 
relative importance of the variables in terms of the desired 
response. Since u and y are the deviations in control and 
output vectors, the minimization of the performance index 
results in minimizing the control effort as v/ell as the 
error in the output response of the system after an initial 
disturbance . 


The problem formulated above does not include any 
constraints of either the equality or the inequality type. 
Usually depending upon the physical nature of the system, 
certain variables are constrained. Optimal regulator 
problems with constraints are difficult to solve. The 
constz'aints can be indirectly taken into account by 
assigning suitable penalties to the constrained variables. 
The Q and R matrices themselves can be chosen to reflect 
these penalties . 

Once the choice of the performance index is made, 
the optimal control law can be determinoJ by the application 
of Pontryagin's maximum principle or the Hamilton-Jacobi 
theory. The optimal control u minimizes the objective 
functional given by the equation (4.2.1), and subject to 
the differential equality constraint 
X = A X + B u 


( 4 . 2 . 2 ) 
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v'/liich rtpresents the state iquatjon of the system and the 
algebraio equality constraint 

Z = 1 '£ ( 4 -. 2 -. 3 ) 

\vhich is the output equation of the system. Where A is a 
nxn matrix, B is a nxr matrix and £ is a pxn matrix < 

The existence of this requires prior investiga- 
tion before a formal solution to the above problem is 
attempted. It is ascertained by invoking the notion of 
output controllability. In the problem under consideration 
.it is desired to transfer the system to its normal op.ersting 
condition (origin fora linearized model) after an initial 
disturbance. The condition of controllability ensures that 
such a transfer of liie system can be affected,. , The presence 
of unGont3;‘,ol Table dynamic inodes implies that for this system, 
bhcre does non exist any control which provides a transfer 
to -^he desired operating condition. 

A system is said to be completely output contro.lla- 
ble if it is possible to construct a control vector 
which will transfer any given initial output ); to any 

final desireci output ,y;(t.|) in, a finite tim C: interval 
t,, ;g: t '^.ti , The condition for the system to be output 
controllable is that the matrix 

GAB OA^B ... G A^”'' .B) 

be of rank p. 



To check the controllability of the system under 
consideration 5 the 3x|l|. matrix P is formed, and the 
determinants of different combinations of 3 columns of 
matrix P were found to be non-zero. Thus the rank of 
P is 3 thereby implying that the system is controllable. 

As a further verification of this result, the system has 
also been reduced to its normal form, and it was found that 
the input matrix has no rows which have all zero elements. 

17 

This ensures that the system is indeed output controllable 

4 . 3 Deriv ation of the Optimal Control Law ; 

Since the theory of optimal output regulator is 
v^ell kno’vn, only the major steps are presented below, for 
the system given by equation (4.2.2) and (4.2.3) and the 
performance index given by equation (4.2.1). 

The Hamiltonian H is first defined as 

H = ^ Q + -I u^ R u + (A X + B u) (4.3.1) 

where A is the Lagrange multiplier vector. 

Por an optimal control u*, X and x are the solutions of 
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that is, 

X = A X + B u 
\ = -C^ i C S - 2 

and S u + B^X = 0 (4. 3-. 3) 

Solving for u, 

u = -R""^ X (4.3.4) 

Assuming that the lagrange multiplier 2 ^ related by 

the matrix equation 

X = K X (4.3.5) 

the optimajk control becomes 

u* = -R"'' K X (4.3.6) 

In the infinite-time case, the symmetric matrix K 
is the solution of the nonlinear Riccati equation 

-i^K-KA + KBE"”' B^K-C^QC = 0 (4.3-7) 

The optimally controlled system dynamics are governed by the 
homogeneous state equation *. 

2 = (4 - B R~'' B^ K) X (4.3.8) 

Brora the above equation it is obvious that the system 
dynamics are altered, hence even an unstable system can be 
stabilized with the application of the optimal control. 
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4 .4 Applioation of Successive Approximation for the 
Solution of the Ricoati Equation ; 

ThG solution of the 7x7 matrix Eiccati equation 

implies the solving of 28 simultaneous nonlinear equations. 

Ihey are difficult to solve analytically and therefore the 

method of Successive Approximation -.as given by Puri and 
16 

Gruver is adopted for the solution. In this method the 
approximation in control policy space is combined with 
stability considerations from the Pirect Method of Liapunov, 
to genexate a sequence of suboptimal control functions, which 
have a monotonic convergence. A ccntrol function is first 
chosen which represents a . stable system. The Hamiltonian 
is then minimized to obtain the second approximation. 

The Ricoati equation at the kth iteration is given by 


= 0 


vdiere 


= A - B R-'' ^ 

and 

qU) 

= Q C + ^ B R' 


The above equation can be solved in a vector equation form 

(4.4.2) 

where = (k^^^ ... k^^^ k^|^ ... k^^^ ... 

S- Cff - iLl J. vectors 



36 




q(k) q(k) 

^1n ^22 


) ^ELl^lUyec tors 


and are the elements of and 

-L J X J 


and is an 


X 5 i.. ^ ± l -l matrix whose elements are the 


GOn'Dination of the a^ ^ of the matrix V . The elements of 

-i* J 

r ic ') 

are obtained by equating the coefficients of 
equation (4.4.1) and is reduced to the followirg form 


,(k) j^cic: 


(4.4.3) 


The solution of vector results in the gain matrix 


The weighting matrices Q and R on the output and 
input vectors are chosen as Q = diag(100 100 100), 

H = diag(1 1). As the technique requires to make the 

closed loop system stable, a choice of K^°^is made and the 
new gain matrix E is then determined. The method of 
Successive Approximation is ccntinued until t he difference 
in the coefficients of matrix K between two consecutive 
iteratior^is less than Mo. The method converges in 10 
iterations for the operating conditions given in Section 2.4. 
The matrices B>A>0 for this operating conditions are 
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- 2.1255 

0.6642 

1.3285 

0 

0 

0 

0 

36.231 

- 76.085 

56.231 

583.14 

- 201 .86 

107.25 

- 0.6540 

24.154 

12.077 

- 38.646 

0 

0 

0 

0 

181.15 

- 694.42 

181 .15 

- 53.829 

40.371 

- 295.12 

- 0.5406 

0 

0 

0 

26 . 9 H 

- 35.886 

0 

0 

0 

0 

0 

0 

0 

0 

1 

- 0.2096 

L 

O . 3 o 39 

- 0.2096 

0.3212 

- 0.1931 

0 

- 0.1667 

i 0.4460 

- 0,9366 

0.4460 

0.5751 

- 0.4313 

- 0.3389 

, 0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

1 


and the natrix K found by the above technique is given by 


5.2204 

- 0.1364 

0.3866 

- 0.0125 

- 0.2300 

2.7671 

- 4.5565 

- 0.1364 

0.4925 

- 0.1455 

0.0065 

0.0013 

0.7798 

- 0.0298 

0.3866 

- 0.1458 

0.0953 

- 0.0088 

- 0.0795 

- 0.0124 

- 0.4965 

- 0.0125 

0.0065 

- 0.0088 

0.3347 

- 0,0866 

0.0762 

0.0328 

- 0.2301 

0.0015 

- 0.0795 

- 0.0866 

0.1565 

- 3.7236 

- 0.4176 

2.7671 

0.7798 

- 0.0124 

0.0762 

- 3.7236 

269.17 

91.553 

-4 - 3565 

- 0.0298 

- 0.4965 

0 . 0328 , 

- 0.4176 

91 .553 

95.133 


4 . 5 Effect of Different Performance Indices and Operating 
Points : 

Different performance indices are generated by keeping 
the R matrix fixed as R = diag(1 1), and the weightages on the 
output vector sequentially being increased, the sequence being 
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Q = ctiag(1 1 1), Q = diag (10 10 10) and Q = diag(100 100 100). 
The output response for these three cases is obtained by 
solving equations(4 .3.8) and (4.2.3) for an initial perturba- 
tion = 0.05 p.u. is shovm in Pigs. 4.1 - 4.3, the 

optimum control u* being aetermined from equation (4.3*6). It 
indicates a progressing increase in the control effort and the 
outputs have faster exponential decays in the increasing 
sequenceP'or the last case the system response dies down within 
2.5 seconds, although at the expense of extra control effort. 
The optimal control at two other operating conditions given 
in Table below 


; Case (ii) 

Case (iii) 

r 

0 1 

O 

O 

10° 

T? 

^fd 

1 .0 

1.0 

V 

1 ,0 

1 .0 

0 

0.9671 

0.9956 

} 

^qo 

-0.2496 

-0.1264 

«t 

0.9990 ' 

1.0710 


is determined. The response of the regulated system for 
an initial perturbation inA"^^ = 0.05 p.u. is shown in 
Pigs- 4.4 and 4-5 reveals that the optimum control effort 
required is less as the load angle decreases. This ought to 
be expected since the restoring torque is more at smaller l...ad 
angles. 
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4 . 6 Validity of the Linear Model ; 

The optimum control law obtained in equation (4.3-6) 
is for the linccrized systera model around a specific operating 
condition. The actual representation of the system model is 
nonlinear as indicated in Chapter II. Thus the validity of 
the optimum control derived on basis of the linearized model 
when applied to the actual nonlinear system remains to be 
investigated. The nonlinear model is represented by 
equations (2,2.29) to (2.2.38) to which the optimum controllaw 
given by the equation (4.3.6) is applied. As the optimum 
control law is represented in terras of the perturbed state 
vector, the latter is obtained by subtracting the actual 
value from the one corresponding to the steady state 
operatirg condition. The resultant output of the controller 
is added to the nominal c cntrol inputs corresponding to tie 
steady state operating condition. This provides the total 
control to the nonlinear system. 

The response for the operating condition given in 
Section 2.4 for the nonlinear model with the optimal control 
law was obtained for an initial perturbation in '^^=0.05 p.u. 
is shown in Jig. 4. 6. The comparison of the response with 
the response obtained for the linearized model indicates, 
that the optimum control determined on the basis of the 
linearized model when implemented on the actual nonlinear 



system gives quite satisfactory performance for small 
perturbation. 

4.7 Gonclusion ; 

The comparison of the response of the system with 
a conventional regulator (Pig. 3. 3 to an optimal regulator 
Pig. 4. 3) for the same perturbation in the variable 
reveals the superiority of the latter. The output response 
deoa/s exponentially, without any overshoot in the optimal 
case . 


In the design of the Optimal Regulator the 

O 

configuration of the controller is not assumed Yu etal 
assume the configuration and have obtained a specific 
optimum control lav. This is an unnecessary constraint on 
their design procedure. 

The optimal control law obtained in equation 
(4. 3* 6) requires the entire state vector for the purpose 
of feedback. The states by themselves are not unique and 
depe-nd upon the particular formulation of the state space 
model. Quite often these state variables are not available 
for direct measurement. On the other hand the system output 
are the ones which are available for control purposes. 
Therefore it becomes necessary that only outputs which are 
directly available should be used in the implementation of 
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tile conUroL law. This is carried out by reconstructing 

. 'K, VA I ‘u, ■ • . „ 

the state Vector frora the available system inputs and outputs 
and using the reconstructed state vector for the control 
purpose 
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CHAPTER Y 

OBSERVER RESIGN FOR THE OPTIMAL 
OUTHJT REGUIA.TOR 


5 . 1 Introduction ; 

The design of the cptimum output regulator in the 
previous chapter assumed the availability of all the state 
variables for the purpose of feedback. It is seldom that 
all of these variables are accessible for direct measure- 
ment. To obviate this difficulty a dynamic compatible 
observer is designed which reconstructs these inaccessible 
state variables from the system input and output variables. 

The reconstructed variables are subject to an exponentially 

1 8 

decaying error. Bongiorno and Youla have given a design 
procedure for the construction of a compatible observer. 

The Optimum Control law derived in e<3.uation (4.3-6) 
Y^hen cascaded with the observer gives the overall configur- 
ation of the optimum regulator. 

5.2 Eormulation of- the Observer Problem ; 

The observer problem is the reconstruction of the 
state vector for the system given by equations (4.2.2) and 
(4.2.3) from the input and output vectors of the system. 
Before the formal reconstruction is considered, it is 
worthwhile to investigate whether the system dynamic modes 
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of behaviour as represented by the state vector can be 
ascertained from measurements of the available inputs and 
outputs. The concept of observability is useful for this 
purpose . 

A system is said to be completely observable in 
the interval t^^ "t^t^ for every t^ and some t.j , if every 
state x('t) can be determined from the knowledge of 
y_(t) and u(t) ±n t^^ t-5^t^. The condition for the system 
to be completely observable is that the composite matrix 

Z = (Ah“-’ ch 

is of rank n. 

To check the observability of the Power System 
under consideration, the 7x21 matrix P is formed for the 
operating conditions described in Section 2.4, and the 
determ j.nants of different c ombin afions of 7 columns of 
matrix P were found to be nonzero. Thus the rank of 
P is 7, thereby implying that the system is observable. 

A.s a further verification of the result, the syston has also 
been reduced to its normal form, and it was found that 
the output matrix has no columns which have all zero 
elements. This ensures that the system is indeed 
observable"''^. 
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Thus for the linear observable system described 
by equations (4.2.2) and (4.2.3), the dynamic observer is 
defined as 

z = llz4-Gy;+Hu 

and £ = ii1 21 + Ii2 ^ (5.2.1) 

where F, G, H, L-] and I2 are real constant (n-p)x(n-p), ,■ 
(n-p) X p, (n--p)x r, nxp nx(n-p) tnctriccs respectively, 

A. 

X is the estimated value of x and 

S. ~ (5.2.2) 

where T is the unique solution of 

TA-FT=G0 (5.2.3) 

provided A and F do not have common eigenvalues. For a 
choice of 

H = i I (5.2.4) 

The error between the actual and observed state vectors at 
time t is given by 

£(t) = exp(F(t - t^)) l(tQ) (5.2.5) 

and the reconstructed state vector is given by 

X = 1(W X + (5.2.6) 

where the vector £ is defined by 


”0 


0 

-p 


-P 

e 


exp(F(t - t^) 0(1^)) 


e 


(5.2.7) 
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In the above equation, 0 ^ denotes the p dimensional column 
'■/ector whose eleoicnts are zero, F is a stable matrix and L 
and t/ ;~re respectively, the partitioned matrices 


ill —2 


( 5 . 2 . 8 ) 


Vv h ere 


- 00 ^( 1 ) 1 

a y 


( 5 . 2 . 9 ) 

( 5 . 2 . 10 ) 


(G F G F G 


“ • if'"'* 2) 






T' OC, withOC = 1 


is the characteristic equation of the matrix A andOC^’s 
are the co-efficients of the characxeristic equation, and 


o 6 i 

ip 

^2 ip 

^n -1 

ip 

/-'f 

C ^2 

ip 

0^3 ip 

OC n 

ip 

* 

cc 

n- 

lip 

• 

o 

OCn ip 

9 tt « % 

0 

• 

% 

OCn 

ip 

% 

* « • • 

2 p 


where I_p and Op are pxp identify and null metrices 
respectively. 



An observer is said to be compatible if its output 
elusls the state of the plant, to within an exponentially 
decaying error. A necessary and sufficient condition for an 
observer to be compatible is the nonsingularity of the 
matrix W given in (5.2.9 )• 

5 . 3 Gonstruction of an Observer ; 

The system under consideration requires the 

c ons; truction of seven variables from three outputs and two 

procedure 

inputs of the plant. Prom the/outlined in Section 5-2, the 
dynamics of the observer should be of the fourth order. The 
matrices P and G of the compatible observer are so chosen 
that the pair (P, G) represents a controllable system and 
the matrix E should not have any eigenvalues in common with 
the original system. Prom the above considerations the F 
and G matrices are chosen as follows: 




— 


— 


— 

o 

^ — 

C\i 

1 

0 

0 


1 

0 

0 

0 -7.385 

0 

0 

G = : 

0 

1 

0 

0 0 

-7.912 

0 


0 

0 

1 

0 0 

0 

-7.342 


0 

0 

1 


The matrices T and W are determined from equations (5.2.10) 

and (5.2.S) respectively and the matrix W is then checked 

for nonsingularity. The matrices H and h are obtained from 

-1 

the relations H = T B and 1 = W . The computation of 



53 


raatrices H, T, and L 2 are done in Double Precision to 
minimize truncation errors, as the parameters of the 
observer are highl-' sensitive, and slight discrepancy leads 
to improper reconstruction of states. The matrices obtained 
are 




0 . 01 9505 

0.000048 





-0.000252 

- 0 . 0031 16 



H 








0.001623 

0.021459 





0.001886 

0.023152 



.56x1 0 "' 

.46x10^ 

-.38x10® 

.13x10® .14xlS^ 

-.10x10® .29x10^ 

-.73x10^ 

.64x10^ 

-4 

.52x10^ - 

.86x10^ .69x10^ 

.13497 -.19x10® 


.47x10^ -.44x10^ -.38x10^ -.49x10^ .28x15^ .12875 

• 55x10'^ -.48x10^ -.38x10^ -.64x1'D^ -.51x1^^ .33xlD^ .13891 




-0.8295x10""' 

0.8454x10 

0.1226x10 

0.1 394x10^ 
-0.9027963 
o' 


0.3568x10^ -0.1287x10^ 

-0.4529x10^ 0.1427x10^ 

-0.5784x10® 0.1927x10^ 

-0.6608x10^ 0.2083x10^ 

0.4625x10® -0.1357x10'^ 
1 0 


0 


0 


1 
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-0. 1431x10^ 

-0.2635x1 0^ 

-0.3205x10^ 

-0.3249x10^ 

0.1846x10"^ 

0.3345x10'“^ 

0.4041x10® 

0.4127x10^ 

0.1953x10^ 

0.4272x10^ 

0.5144x10^^ 

0.5273x10® 

0.2765x10^ 

0.4880x1 o’* ° 

0.5897x10® 

0.6021x10^ 

-0.1341x10^ 

-0,3415x10^ 

- 0 . 4122 x 1 o’^ 

-0.4215x10® 

0 

0 

0 

0 

0 

0 

0 

0 


. 4 Perforaance of the Optimum Output Feedback Regulator ; 

1 8 

Borgi^omo and Youla have pointed out .that there 
is an increase in the values of the performance index when 
a compatible observer is employed to reconstruct the state 
vector, which in turn is employed as the input to .the optimal 
feedback regulator. Furthermore they have shown that there 
exists no other choice of the feedback gain matrix which will 
give a smaller increase in the value of the performance index 
for all possible initial conditions, when the practical 
requirement of stability is imposed. Thus the observer 
designed in Section 5.3 can be introduced into the ^stem 
which enables the implementation of the optimal control with 
output feedback. The control is then fed to the unregulated 
system. 
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The output response of the system for different 
initial perturbations of 0.05 p.u. in each of , 

Z\y''kd» A'y'q* ® obtained. 

Figures (5.1 to 5.7 show the resulting performance curves, 
for iiie case when the initial condition vector for the 
observer is exact. The error introduced in the overall 
response due to the introduction of the observer is 
negligible .and this is seen by comparing Figure 4.3 with 
Figure 5.1 for the same initial perturbation p.u. 


The response due to initial perturbations in^"^^, 
andA^^q do not have any marked effect on the system 


performance. This is due to the fact that the time constant 


associated with the damper coil in the quadrature axis is 
small. Perturbation in ^5^ andA52 requires a large 
controlling effort in A T^ as compared to other perturba- 
tions, to brirg the system to the steady state operating 


condition quickly. This may be attributed to the inertia 


of the generator which is to be overcome after the initial 
perturbation. The perturbation requires a greater 
change in field excitation control rather than the input 


torquta, as the field excitation has a major effect on the 


flux linkage in the field winding. In the case of 
perturbation inA'^^ the control effort from both the 
inputs are almost equal. In all these cases the output 
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response dies down exponentially and the steady state 
operating point is reached fairly fast within 2.5 to 3*0 secs, 
of the initial perturl^ation in the state variables. 

The output response with zero initial conditions 
for the observer states, for an initial perturbation in 

d “ ® .05 p.u. is shown in Pig. 5. 8. It differs from the 
output response y/ith exact initial conditions for the observer 
during the initial stages of the response and takes a slightly 
longer time to reach the steady state operating condition. 
Thus the improper choice of the initial conditions for the 
observer states has an adverse effect during the initial 
stages of the response only. 

The differences in the control law expressions 
obtained in Section 4.6 for three different operating 
conditions are relatively small. Thus the same optimal 
output regulator configuration can be used for all these 
three operating conditions. However, the same thing cannot 
be done regarding the observer. The observer is highly 
sensitive with respect to the changes in the system 
co-efficient matrices as well as the initial conditions, and 
therefore the observer designed for one operating condition 
when used for a different operating condition wlill have 
adverse effect on the system response due to the improper 
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reconstruction of the state. Thus care should be taken 
to always use a proper observer for each operating condition 
of the system. 

When the optimal routput regulator designed for 
ease (ii) is used along with the observer designed for each 
specific operating condition, the response for such a 
combination for the other two operating ccnditions are shown 
in Pigs. 5. 9 and 5.10. It is slightly inferior as compared 
to optimal output regulator and observer for the specific 
operating condition, thus justifying the use of a common 
control law. 

5 .5 Transfer hatrix Bepresentstion of the Optimal Output 
Peedback Piegulator ; 

The observer and the optimal control law cascaded 
together can be thought of as constituting the Optimal 
Output Peedback Regulator. The observer discussed in 
Section 5.2 is given by 

z = Pz + Gy;4-Hu 

and % = hi IL 112 ^ (5.4.1) 

the optimum control law as obtained in Section 4.2 is given 

by 

(5. 4. 2) 

-1 T 

P = -R B K 


where 
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It is observed in Section 5 .3 that the effect of assuming 
zero initial conditions for the observer alters only the 
initial portion of the response, the error due to improper 
initial conditions is reasonably small. In view of this 
it is convenient to assume the initial conditions to be zero. 
Using this fact, and taking the Laplace Transform of the 
eluotions (5.4.1) and (5.4.2) and rearranging 

(s I - F) Z(s) = G Y(s) + H U(s) 
or Z(s) = (s I - (G Y(s) + H U(s) ) 

Z(s) = L(s) Y(s) + M(s) IJ(s) (5.4.3) 

where 

D(s) = (s I - F)"”' G 
and M(s) = (si- F)~^ H 

X(s) = Y(s) + I 2 Z(s) 

and U(s) = P X(s) (5.4.4) 

But X(s) X(s) 

U(s) = P(s) Y(s) + P(s) L 2 L(s) Y(s) + P(s) Ig M(s)U(s) 
or U(s) = U(s) Y(s) (5.4.5) 

where 

]i(s) = (I - P I 12 M(s))“^ (P + P Ig D(s)) 

In the above equation, the matrices P, , L 2 , M(s) and D(s) 
are all known and U(s) can be calculated. 
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A typical configuration of the Optimal Output 
Eegulator for the operating condition given in Section 2.4 
is shown in Fig.S.II. 'The implementation of the Optimal 
Output Regulator thus reduces to the realization of the six 
elements of the 2x3 transfer matrix F(s), shown in Rig.5.11. 

The outputs of the controller are fed to an exciter to 
develop the required change in field excitation md to a 
motor to modify the governor setting to produce the required 
change in the input torque. Thus, the dynamics of the 
exciter and the motor must be included in realizing elements 
of the transfer matrix ^(s). The rest of the dynamics can 
be synthesized using appropriate compensating devices. The 
inputs to the optimal output regulator can be obtained by 
the following technique. The change in the derivative of 
the load angle can bo obtained as proportional to the 
difference between a fixed reference and the actual tacho- 
generator output and integration of the resultant gives the 
change in load angle. The difference in the terminal voltage 
and the fixed reference voltage gives the other input. 

The outputs of the optimal regulating equipment can 
be made effective as follows. Charge in excitation voltage 
can be implemented by directly changing the excitation voltage. 
The other change in the input torque is achieved by ccn trolling 
the gate opening in the case of a hydro system, or by chang- 
ing the amount of steam to the turbine in a thermal system. 








UC6 p*5 



Pp,,6t3 FfirfoTmajnco of th-e systera with t.b© optimal output foo<iao, 
rerulatoT tm case (1) with exact .initial cor iii 1 oi;:; .fc »■ 
the observer. 
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may include them as they interact with the generator dynamics. 
An optimal regulator can be designed for the overall system 
controlling the prime mover also along v/ith the generator. 

(iii) In this thesis the problem formulated is solved 
as an infinite time output regulator problem. Other possible 
performance objectives such as minimum time with constraints 
on the magnitudes of 'the input variables can also be investi-" 
gated, although these are difficult to solve analytically 
Besides, 1he minimum time criterion usually results in a 
bang-bang control, which cells for switching between the 
extreme values of the permissible control magnitude. Such a 
large variation in control may render the linearized model 
chosen for the genertlor invalid. 

(iv) The theory of observers under certain conditions 
ensures the rec aistruction of the state vector within a small 
error. The design procedure, however, cannot incorporate 
information regarding the practical aspects of the problem. 

The observer should have smaller eigenvalues in order that , 
its hardware realization is feasible but this has the drawback 
that the error in the reconstructed state tends to get large. 

A suitable compr.cmise has to be made between these conflicting 
rcQuirements. 

(v) In certain situations it is advantageous to use 

a simplified form of the optimal control law in actual 'w 

implementation. This simplification can be performed by the 

19 

method of aggregation suggested by Aoki . 
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APPENDIX A 

DESIYATION OE LINEABIZED MODEL 


In the derivation of the linearized model the 
technique of linearization about an operating point is 
applied, by substituting the perturbed values for the 
variable quantities in System equations (2.2.29) to (2.2.38). 

The perturbed equations are; 


PCFf^inW-pJ = 

^f1 


"^o^fd 


do^^Tfd^ - o (pf a fd ta o^fa ) 


’^O^fd 


^fdko^kd'^ x_^ (E^^q+Ae^^^) 




ad 


"fd 


(A.1) 


vv^(r+e) vvQ(r+e) 


P%oWd^ “ K.jX£.,XQ.,C^fdo^i?^fd^ ^ 
w^(r+e) 

, 4.Mr -1— r — I — . 




^l^al^kdl 


^a1 Va1 


-f 


^"2^a1^kq1 


:p^^'^qo+i^kq^+ w^gvSin(A^ q-^S ^ )+w^hvGos(<5 ^ q4^ ^ ) 

..... (A.2) 


P^Kdo^^kd ^ 1 ^fdo^^d ^ 

^kdl *'^"l ^kd1 (-^-‘5) 
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+ _ 2 l-ri_ - 1 - 
^a 1 ' 


Tdo^'fe' 


^l^al^kdl 




o'^?kd 


(r+e) 


w^(r+e) 


^KoX^r”' ^ K 2 Xg ^x^q^ 


"a 1 "^ 2 ^a 1 


Y^^hYSlXl(S^Q+hSl) + 


Gos( 6 ’^q+A<S| ) 


(A. 4 ) 




^ 2 ^kq 1 ^a 1 


^ql "^ 2 ^kqi 


•••*(A* 5 ) 


P "I Q'^Ac.i-j ; — ( ^^*20 ^ 


(A. 6) 


P^^20'^A42) ~ ji ^ '^i”^^d ^^20"^ ^2^' 


i'^A o+ 


^qo^A^K) 




^Fd o+tt ) C%^a o'^^lcd ^ Cfqo^M ^ ^fd o"^ ^ 


"' 2 ^a 1 ^kq 1 






K,Xg,Xji 


^l^al^kdl 


(A. 7) 


(e_ +/de+)^ = (e.^+ie^)^ + (e„^ 


do^^'^d 


^qo 


(A. 8) 


( 6 1 o'^A^i ) = (tj^Q +/iS^ ) 


(A.9) 


( b 20 '*'A^^ 2 ^ ~ 


(A.IO) 
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(A. 11) 


(A. 12) 


Eliminating the steady state terms defined in e*luations 
(2.2.29) to (2.2.38) and retaining first order terms results 
in the foilov‘/ing linearized equations relating the perturbed 
quantities. 


Hftd 


v/„r 


p-"fd, 1 




-.1 


w_r , 


fd Kix^.x„^ vd^ K^x, . .xvi-^rkd 


'-l^fl^al 


^l^kdnfi' 




(A. 13) 


WQ(r+e) 




-f 


w^{r+e) , 


X 


a1 


^1^1 


-1 )d^' 


'r 

d 


AO/ , '^‘V 

^kd + 


^1^a1^kd1 


1 w f 


O X 


a1 




+ (w^gvOos) ho “ w^hvSin5^Q)/ld^ ^T’qc^^^a 


(A.14) 
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kd 


0 ka 


l\(t 


^o^'kd 


, %^kd/ 1 

- + — r U 


*. ..(A/15) 


^A% = 


w^f ^ , w^f ^ 

•4 fed - + -f-T^ 




° ^a1 ^l^al 


-1 




’1""a1"^kd1 


A% 


vv'^(r+e) 


w. 


fr+e) 




+ 


2"a1""kq1 


0:. 


( w^hvO 0 Q-^ w^gYSin S^Q)Ad^ -'¥d^^2 


(A. 16) 


A^ '^o"kq 

P/J4^kq = 


Vr 


"'o’^k?, 1 


Vkq1=‘ar'"1 


^kql K2^k41 




kq 


(A. 17) 


'0 /\ /■ .' -1 =4 C 


Paa 1 


(A. 18) 


^dOz 


A% 


^qo kqo __ T^<Lo 


r^f r ^ Va1 ^2^a1 ^^2^kq1 ^ ^l^al 


l/'qo 


Vl^^al^kdl 


■ 


1 / Tdo ^ "^do Tfdo yi&o 


^E^al ^2^a1 ^1^a1 


Kl^f1 ^ ^l^kdl'^^^ 




T 


(A. 19) 


D 


Substituting (A. 11 ) 9 nd(A^ 1 ' 2 ) in (A^ 8 ), and simplifying 




= o^'-nf1 a1 
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+ 




-- ( "I 


®oKl^a1^kd1 ®o^a1 Va7 






“ ^sind 


(A. 20) 


Adj 

=~as ^ 

(A. 21) 

A52 

=Zi6^2 

(A. 22) 



APPEEDIX B 


RUHGB igJTTA METHOI) FOR S GLUT I OF OF 
DI EPEREMTIAL EQUATI OPT S 


Range Kutta fourth order method of solvi% the 
given system of differential equations is described as 
f ollovvs; 


Eor the system (X) = f(X, t) 


(X) at t + A't can be found from (X) at t as 


t+At — t 


+-2X2 + 2 K^ + 
5 


where 


= 

t f(x,t) 


K2 = 

■fc ICX + t , 

tv+ 

II 

t f(x + i K2, 

t + %) 

II 

t :f(X + K3 , 

t + 4 t') 


(B.l) 


The above equation is obtained by using the truncated 
Taylor series expansion of the function about the initial 

point X^ and equating the corresponding coefficients- of 

n 

(/^t) , where n = 1,2,3»4. This method has rounding of 

5 

errors proportional to (At) and is numerically stable for 
reasonable values of At. 



APPsiroix c 


The value of the rotor inertia time constant (T ) 

m 

of the generator (p.16) is somewhat large. Instead if the 

value of the inertia constant is tekan as 0.019 seconds as 

14 

taken by Prabhashankar etal , while the other data of the 
system given in Section 2.4 remains unaltered, the following 
are the corresponding changes: 

In equation 2.3-1 on p.36 


B = 


0.3454 

0 


0 

0 


0 

0 


0 

0 


0 

0 


0 

0 


A _ 


In equation 2.3-1 on p.-37 
- 2,1255 0.6642 1.3285 0 


0 


0 

52.333 

0 


T 


0 


36.231 -76.085 36,231 583.14 -201.86 107.25 -0^6540 


24.,154. 12.077 -38.646 

181.150 -694.42 181 .15 
0 0 0 

0 0 0 

-65.819 95.438 -65.819 


0 


0 


0 


0 


-53.829 40.371 -295.12 -0.5406 

26.914 -35.886 0 0 

0 0 0 1 


100.87 -60.624 


0 -52.333 



In eHuation 4.3.7 on p.37 


5.2642 
-0.1394 
0.3932 
-0.0109 
- 0.221 3 
2.7o39 
- 0,0078 


-0.1394 

0.5027 

-0,1467 

0.0050 

0.0025 

0.6941 

-0,0025 


0.3932 

-0.1467 

0.1025 

-0,0097 

-0.0742 

0.2037 

-0.0114 


-0.0109 

0.0050 

-0.0097 

0.3470 

-0.0967 

0.0918 

0.0264 


-0.2213 
0.0025 
-0.0742 
-0.0967 
0.1597 
-2. 7176 
-0.0225 


2.7839 

0.6941 

0.2057 

0.0918 

-2.7176 

152.90 

0.2738 


-0.0078 

-0.0025 

-0.0114 

0.0264 

-0.0225 

0,2738 

0.1734 


in equation 5.2.1 on p. 52 


I — 

-7.50 

0 

0 

0 


5.0 

0 

5.0 r 


-7.50 

0 

0 

G -= 

^ 1 .0 

1 .0 

0 

0 

0 

-7.50 

0 


■ 1 ,0 

1 .0 

0 

4 0 

0 

0 

-7.50 


1 1 -0 

0 

1 .0 


In eqaation 5.24 on p. 53 


H 


-46.171 -6.9687 

0.3339 0.1827 

0.3339 


-0.9234 


0.1827 
-1 .3937 





In equation t.2.10 on p. 53 


-1 .3367 

.17x10‘' 

.98xlS’’ 

-. 13 x 1 ^’’ 

.22411 

-.96604 

-.13316 i: 

0.97x10 

.74x10'^ 

-.71x10^ 

.18x10^ 

.98x1D^ 

.15755 

.35x10^[ 

0.97x10 

.74x15^ 

-.71x10^ 

.18x10^ 

.98x10^ 

.15755 

.35x10^ 

-.26755 

,33x10^ 

. 20x10*^ 

-.26x10^ 

.45x10^ 

-.19321 

-,27x10*' 


In equation 5.2.8 

on p.53 





n 







^ -0.2057 

-0.6036 

-0.1626 





0.2274 

-3.4,560 

0.2037 





-2 . 0465 

7.0544 

-1 .1136 , 




5ii = 

3.7533 

-12.112 

1 .4549 





-0.1 365 

-2.8041 

0.1778 





0 

1 

0 ■ 

I 

[ 




0 

0 

1 




In equation 5. 

2 on p.54 






r~ 


-0.6567x10*'^ 

0.1989x10**^ 

-0.1989x10*''^ 

0.3284x10*'® 

0.6334x10*”^ 

-0.5797x1 o’’ 

0.5797x10*''^ 

-0.3167x10*'® 

-0.7526x10*'^ 

0.2190x10*'® 

-0.2190x10*'® 

0.3765x1 0*'’^ 

0.7174x10*'^ 

-0.2708x10*'® 

0.2708x10*'® 

-0.3587x10*''^ 

-0.2706x10*'^ 

0.1189x10^^ 

-0.1189x10^"^ 

0. 1353x10*'® 

0 

0 

0 

0 


0 


0 


0 


0 




The transfer matrix N(s) defined by equation 5.4.5 on p.71 is 







When the exact initial conditions are not known and 
therefore soiDe arbitrary values are assumed, Ihe choice of the 
observer dynamics play a key role in determining how fast the 
initial condition error response decays. For this reason, it is 
desirable to keep the eigenvalues of F matrix reasonably large. 
On the otherhand , choice of the observer matrices F and G both 
play an important role in theoverall configuration of the 
optimal output regulator. ?/ith different choices these matrices 
so long as the observer has faster dynamics only the initial 
portion of the overall output response of the system gets 
affected, 

The qualitative performance of the system with the 
optimal output feedback r egulator , with the new inertia time 
constant, is essentially the some except for an increase in the 
input torque required during the initial stages of the response. 
The effect on the Tran sfer Matrix Representation of the Optimal 
Output feedback controller is however much more pronounced. The 
D.C. gains associated with the individual transfer functions 
are no?v within a reasonable range as compared to ihe large 
values of these gains earlier. 
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